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ABSTRACT 

We follow the chemical evolution of the Galaxy for elements from Ba to Eu, 
using an evolutionary model suitable to reproduce a large set of Galactic (local 
and non local) and extragalactic constraints. Input stellar yields for neutron-rich 
nuclei have been separated into their s-process and r-process components. The 
production of s-process elements in thermally pulsing asymptotic giant branch 
stars of low mass proceeds from the combined operation of two neutron sources: 
the dominant reaction 13 C(o;,n) 16 0, which releases neutrons in radiative condi- 
tions during the interpulse phase, and the reaction 22 Ne(a,n) 25 Mg, marginally 
activated during thermal instabilities. The resulting s-process distribution is 
strongly dependent on the stellar metallicity. For the standard model discussed 
in this paper, it shows a sharp production of the Ba-peak elements around 
Z ~ Zq/A. Concerning the r-process yields, we assume that the production 
of r-nuclei is a primary process occurring in stars near the lowest mass limit for 
Type II supernova progenitors. The r-contribution to each nucleus is computed 
as the difference between its solar abundance and its s-contribution given by the 
Galactic chemical evolution model at the epoch of the solar system formation. 
We compare our results with spectroscopic abundances of elements from Ba to 
Eu at various metallicities (mainly from F and G stars) showing that the ob- 
served trends can be understood in the light of the present knowledge of neutron 
capture nucleosynthesis. Finally, we discuss a number of emerging features that 
deserve further scrutiny. 

Subject headings: nucleosynthesis - stars: abundances, AGB and post-AGB - 
Galaxy: evolution, abundances 



- 3- 



1. Introduction 

A quantitative understanding of the Galactic evolution of nuclei heavier than iron has 
been so far a challenging problem. Although some of the basic tools of this task have been 
presented several years ago, both from an observational (Spite & Spite 1978, 1979; Sneden 
& Parthasaraty 1983; Gilroy et al. 1988) and a theoretical point of view (Truran 1981), only 
recently the observational data have grown sufficiently in number and precision to allow a 
direct comparison with theoretical predictions (see Mc William 1995 for a review). In order 
to reconstruct the evolutionary history of neutron-rich elements in the Galaxy one has to 
disentangle the s- and r-contributions of all their isotopes and follow their abundance as a 
function of Galactic age. The results must then be compared with the available spectroscopic 
observations of stars at different metallicities, both in the Galactic halo and disk, mainly 
from F and G dwarfs and giants whose surface abundances reflect the composition of the 
gas from which they formed (Gratton & Sneden 1994; Edvardsson et al. 1993). 

Already in the seminal work by Burbidge et al. (1957), the synthesis of neutron-rich 
nuclei was attributed to two different mechanisms. The first (the s-process) is characterized 
by slow neutron captures and occurs mainly during hydrostatic He-burning phases of stellar 
evolution. Here slow means that most unstable nuclei encountered along the neutron capture 
path have time to decay before capturing another neutron. The second mechanism (the re- 
process) is due to the more rapid neutron captures and is generally associated with explosive 
conditions in supernovae (hereafter SNe). For a comprehensive review on the astrophysical 
sites of neutron capture nucleosynthesis, see Wheeler, Sneden, & Truran (1989). 

The current understanding of the s-process is supported by many observational and 
theoretical works (summarized in Sect. 2), indicating that the heavier s-nuclei, from Sr to 
Pb, belonging to the so-called main component (Clayton & Ward 1974; Kappeler et al. 1982), 
are synthesized during the thermally pulsing asymptotic giant branch (TP-AGB) phase of 
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low-mass stars, mainly in the mass range ~ 1-4 M (see e.g. Gallino, Busso, & Lugaro 1997; 
Busso, Gallino, & Wasserburg 1999, hereafter BGW). A different site for the production of 
s-process nuclei is required to account for the weak component, i.e. for s-nuclei below the 
Sr-peak. This site was identified in the advanced evolutionary phases of massive stars (Lamb 
et al. 1977; Prantzos et al. 1990; Raiteri et al. 1993). 

As for the r-process, the actual astrophysical environment is still a matter of debate 
(Woosley et al. 1994; Wheeler, Cowan, & Hillebrandt 1998). It has been recently shown that 
at least two different supernova sources are required for the synthesis of the r-nuclei below 
or beyond the neutron magic number N = 82 (Wasserburg, Busso, & Gallino 1996; see also 
Sneden et al. 1998). For the elements considered here, we make the simplifying hypothesis 
that their r-contributions in the solar system can be considered of primary origin (i.e. with 
yields independent of the initial stellar metallicity). Hence, after the s-contributions at 
the epoch of the solar formation have been estimated with the use of a reliable model for 
the chemical evolution of the Galaxy, the r-process abundance fractions can be derived by 
subtracting the s-fractions from the solar abundances. 

The study of the Galactic evolution of elements produced by s- and r-processes was 
addressed in previous investigations only at a semi-quantitative level. For example, Andreani, 
Vangioni-Flam, & Audouze (1988), and, more recently, Pagel & Tautvaisiene (1997) have 
assumed ad hoc stellar yields and time delays for representative s and r elements in order 
to fit the observational data. With a different approach, Mathews, Bazan, & Cowan (1992), 
following a preliminary analysis by Clayton (1988), examined the effects introduced by the 
primary nature of the 13 C neutron source on the s-process yields, and inferred that for 
metallicities Z > 10~ 3 Z Q the neutron exposure r = J n n VTdt (where n n is the neutron 
density and vt the thermal velocity) roughly scales as Z~ l . Under the crude assumption 
that the production (in mass fraction) of a given s-isotope depends on the product of the 



- 5 - 



neutron exposure times the abundance of iron group seeds, Mathews et al. (1992) estimated 
the stellar s yields by simply scaling the solar s-abundances with metallicity. Today the 
work of Mathews et al. (1992) needs to be revised, as quantitative yields are now available 
from detailed nucleosynthesis calculations based on TP-AGB models for different masses 
and metallicities, with updated reaction networks and neutron capture cross sections. In 
addition, as anticipated by Truran (1981), at low metallicities the r-process contribution to 
the various elements becomes dominant. 

In this paper, we concentrate on elements from Ba to Eu, not only because they belong 
to a major abundance peak, but also because they contain species of very different origin, 
ranging from almost pure r-process production, as is the case for Eu, to a dominant s-process 
origin, as is the case for Ba. The paper is organized as follows: in Sect. 2 we analyze the 
neutron capture process in AGB stars; in Sect. 3 we describe the adopted model for the 
chemical evolution of the Galaxy and the updates introduced in the present work; in Sect. 4 
we show the model results for the elements considered (Ba, La, Ce, Pr, Nd, Sm, and Eu) 
and for their isotopes; finally, in Sect. 5 we summarize the main conclusions and point out 
a few aspects deserving further analysis. 

2. The production of s-process nuclei in Asymptotic Giant Branch stars 

2.1. Previous work 

The first quantitative model for the production of s-process nuclei with A > 85 was 
presented by Truran & Iben (1977) and by Iben & Truran (1978). It was based on the 
activation of the 22 Ne(a,n) 25 Mg neutron source in TP-AGB stars of intermediate mass 
(5 <^ M/Mq <^ 8). However, subsequent reanalyses showed that the attempts to reproduce 
the solar main component through the 22 Ne(a,n) 25 Mg neutron source suffer from major nu- 
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clear problems (Despain 1980; Howard et al. 1986; Mathews et al. 1986; Busso et al. 1988). 
On the other hand, observations revealed that the luminosities of MS, S and C (N type) 
stars are generally much lower than predicted by the above models, thus pointing out that 
the dominant contributors to the synthesis of carbon and s-process elements are AGB stars 
of lower mass (see e.g. Clegg, Lambert, & Bell 1979; Blanco, McCarthy, & Blanco 1980; 
Smith & Lambert 1986; Lambert et al. 1991). In low-mass AGB stars the maximum tem- 
perature reached during thermal instabilties in the He shell is too low (T max ~ 3 x 10 8 K) 
for 22 Ne to be significantly consumed. Consequently, neutron captures must be driven by 
the alternative 13 C(a,n) 16 reaction (Cameron 1954), which can easily be activated at much 
lower temperatures (for the rate of this reaction see Caughlan & Fowler 1988 and Denker et 
al. 1995). 

The amount of 13 C left in the He intershell from previous H burning is by far too low to 
be of any relevance in the neutron production. A challenging problem is therefore to explain 
how a sufficient amount of 13 C can be built locally. The common hypothesis is that some 
protons can penetrate below the formal convective border of the envelope into the top layers 
of the He intershell. At H re-ignition, these protons are captured by the abundant 12 C left 
by partial He burning in the previous convective pulse, giving rise to a tiny layer enriched 
in 13 C (the so-called 13 C pocket). The first mechanism considered in order to generate a 13 C 
pocket was the semi-convective mixing of H into the 12 C-rich and He-rich intershell, driven 
by increased carbon opacity in the cool expanding envelope after the quenching of a thermal 
pulse (Iben & Renzini 1982a,b). This process was shown to be relevant for Population II stars 
(Hollowell & Iben 1988, 1989). Once formed, the 13 C pocket was assumed to be engulfed 
by the next convective pulse releasing neutrons. The resulting neutron density, originally 
estimated to be very high (> 10 11 n/cm 3 , see Malaney 1986) was subsequently shown to be 
of the order of (5-10) xlO 8 n/cm 3 (Kappeler et al 1990; Gallino, Raiteri, & Busso 1993). 
However, semiconvection was found to become inefficient in Galactic disk stars (Iben 1983), 
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and no obvious alternative from simple diffusive mixing was found to work (Iben 1982). 

The penetration of the convective envelope into the partially He-burnt zone after each 
pulse (the so-called third dredge-up) leaves a sharp H/C discontinuity between the H-rich 
envelope and the 12 C-rich and He-rich intershell, which is likely to be smoothed by some sort 
of mixing at the interface, by chemical diffusion, semiconvection, or rotational shear (Langer 
et al. 1999). To address this point, a rather sophisticated treatment of the hydrodynamic 
behavior at the H/He boundary would be required. First attempts with a diffusive approach 
(Herwig et al. 1997), or with a fully hydrodynamical formulation (Singh, Roxburgh, & Chan 
1998) have been recently presented and do show formation of a 13 C reservoir. Indeed, 
penetration of protons at the top of the 12 C-rich intershell appears to be plausible, but 
the mass involved and the resulting H profile have still to be treated as relatively free 
parameters. A dedicated effort on hydrodynamical grounds constitutes a major challenge 
for AGB calculations and s-process nucleosynthesis studies. 



2.2. s-process yields from AGB stars: FRANEC results 

A more indirect approach is to consider different choices for the mass and profiles of the 
13 C pocket, and to estimate their effects on the nucleosynthesis results for AGB models of 
various metallicities. The TP-AGB models by Straniero et al. (1995, 1997) showed that any 
amount of 13 C produced by a proton penetration below the convective envelope border burns 
radiatively in the interpulse period, at a relatively low temperature (T ^ 0.9 x 10 s K) and 
neutron density (n n <^ 10 n/cm 3 ), giving rise to a tiny layer extremely enriched in s-process 
isotopes that will be engulfed by the next convective instability. The resulting pattern of 
AGB nucleosynthesis, and its dependence on the initial metallicity of the star, have been 
recently discussed by Gallino et al. (1999). The changes introduced by these recent studies 
are important: not only the neutron density is lower by a factor of ~ 20 than in previous 
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estimates, but the distribution of neutron exposures can no longer be approximated by an 
exponential distribution, as commonly assumed previously with relevant consequences for 
the dependence of neutron captures on metallicity (Busso et al. 1997). A thorough analysis 
of the implications of the new models is presented by BGW. 

The stellar evolution results shown in this paper are all based on calculations performed 
with FRANEC (Frascati Raphson- Newton Evolutionary Code, see Chieffi & Straniero 1989). 
Stellar models appropriate for Population I stars have already been presented in Straniero 
et al. (1997) and Gallino et al. (1998); additional models have been computed specifically 
for the present paper. Our set of AGB models includes stars of masses from 1.5, 2, 3 and 
5 Mq, and metallicity Z = 0.006, 0.01 and 0.02. Mass-loss was taken into account adopting 
the parametrization by Reimers (1975) with values of the r\ parameter ranging from to 5. 

The FRANEC AGB models reproduce self-consistently the third dredge-up episode and 
the consequent mixing of freshly synthesized s-processed material (together with 4 He and 
12 C) to the surface of the star. According to Straniero et al. (1995, 1997) and Gallino et al. 
(1998), the third dredge-up starts after a limited number of thermal pulses, increasing its 
efficiency up to a maximum while the pulsing mechanism goes on, then decreases as mass 
loss and the advancement of the H-burning shell reduce the envelope mass, and eventually 
ceases when the latter decreases below about 0.5 M & . 

The mass and the profile of the 13 C pocket were constrained by Busso et al. (1995, see 
also Busso et al. 1999) by comparing the results of neutron capture models with spectro- 
scopically determined abundances in Population I AGB stars and in Population II Ba and 
CH stars with [Fe/H] < -1 (Smith & Lambert 1990; Luck & Bond 1991; Vanture 1992; 
Plez, Smith, & Lambert 1993). This was done by keeping constant the total mass of the 
13 C pocket and parametrically varying the concentration of 13 C by factors from to 2 times 
the profile indicated as "standard" (ST) by Gallino et al. (1998), corresponding for the choice 
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of the mass of the 13 C pocket made for this work to 3.1 x 1CT 6 M of 13 C. As a result of this 
study, the ST choice for the 13 C pocket, when used in an AGB stellar model of metallicity 
Z = Zq/2, was shown to reproduce remarkably well the abundance distribution of the solar 
main component (see also Arlandini et al. 1999, in preparation). However, AGB stars in the 
Galactic disk show a significant spread in their abundances of s-process elements (see discus- 
sion in Busso et al. 1995 and Busso et al. 1999). This spread can be attributed, at least in 
part, to variations in the 13 C profile, since the s-process efficiency depends strongly and non 
linearly on the local 13 C abundance. Following this indication, we have explored different 
extensions in mass of the 13 C pocket, while keeping for all metallicities the same 13 C profile 
as in ST case. We will later comment briefly on these choices. Let us mention here that the 
same choice was found to be appropriate to account for the s-process isotopic signatures of 
heavy elements (among which Ba, Nd, and Sm) in presolar grains most likely condensed in 
the circumstellar envelopes of AGB stars of solar metallicity (Gallino et al. 1993; Gallino et 
al. 1997) and recovered in meteoritic material (Ott & Begemann 1990; Prombo et al. 1993; 
Zinner 1997; Hoppe & Ott 1997). 

At the end of the TP- AGB phase, the s-process contributions to the ISM are determined 
by the amount of matter cumulatively dredged-up from the He intershell to the surface and 
lost by stellar winds. This quantity is shown in Tab. 1 for the 2, 3, and 5 M stellar models 
and different metallicities (for the 5 M Q star we did not compute the model at Z = 0.006 
but we doubled the cumulative dredged-up mass mixed into the envelope with respect to the 
solar case, in analogy with the results obtained for the 2 and 3 M Q cases). 

In Fig. 1 we show the abundance by mass of Ba in the He shell material cumulatively 
mixed to the surface as a function of [Fe/H] for a 2 M Q AGB model, for different assumptions 
on the total amount of 13 C nuclei present in the 13 C pocket. The dashed thick lines represent 
the extreme values chosen in the parametrization (0.17 and 1.33 times the ST case), whereas 
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TABLE 1 

Cumulative dredged-up mass in the envelope (in M q ) for AGB stars 



z 


2 Mq 


3 M 


5 Mq 


0.006 


0.030 


0.078 


0.12 


0.02 


0.018 


0.044 


0.06 



thin lines show intermediate cases. The thick continuous line represents the unweighted 
average yield calculated over all the cases shown. In Fig. 2 we show the full set of unweighted 
average yields for Ba, La, Ce, Nd, Pr, Sm, and Eu. 

It is evident from Fig. 1 and 2 that there is a strong dependence of s-process yields 
on the initial stellar metallicity. This important result deserves additional comments. Since 
the build up of heavy nuclei requires neutron captures starting from Fe seeds, the s-process 
is expected to decline with declining metallicity, i.e. to be of secondary nature. However, 
the abundances produced depend not only on the initial Fe concentration, but also on the 
neutron exposure. The concentration of 13 C in the pocket is of primary nature (it is built 
from H and freshly made C, hence is independent of metallicity), while the abundance of 
the neutron absorber 56 Fe varies linearly with Z. Therefore the neutron exposure r (roughly 
proportional to the ratio 13 C/ 56 Fe) is expected to scale roughly as Z~ l . This dependence 
would compensate the mentioned secondary nature of the s-elements, if the yields of s-nuclei 
were linearly dependent on both r and Z. 

One has however to notice that the dependence of s-process yields on the neutron 
exposure is very complex and non linear (Clayton 1968). In general, the production factor of 
a given s-element (e.g. Ba) depends on the average abundances achieved by all the neutron- 
rich nuclei, and varies with Z also for metallicities close to solar. As shown in Fig. 1, in 
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the Galactic disk, for a range of metallicities covering about 1 dex, the strong increase of 
the s-process yields with increasing neutron exposure dominates and the secondary nature 
of these nuclei is overcompensated, so that they actually increase for decreasing Z, instead 
of remaning constant, as the previous simple discussion would have suggested. 

At low metallicities, an important role is also played by light neutron absorbers, which 
are mainly primary in nature, as stressed by Clayton (1988) and by Mathews et al. (1992). 
Among these nuclei one has to include 12 C, 16 and most of 22 Ne with its progeny. This 
is so because a large fraction of 22 Ne ultimately derives from the processing of the primary 
12 C dredged-up into the envelope from the He-shell, which is transformed first to 14 N by 
shell H-burning and then to 18 O and 22 Ne by He-burning. Notice that the recently measured 
neutron capture cross section of one of the major primary neutron absorbers, 16 (Igashira 
et al. 1995), is 170 times higher than the previous theoretical estimate quoted in Beer, Voss, 
& Winters (1992). On the whole, the importance of the moderating effect on the neutron 
density by all light neutron absorbers cannot be discarded. 

In addition, as discussed by Gallino et al. (1998), radiative 13 C burning builds the dis- 
tribution of s-elements in steps of relatively large n-exposition. This gives rise to higher 
values of the Ba/Zr ratio than feasible in any process working inside the pulses (that would 
produce simpler exponential n-exposure distributions). This last feature is a necessary con- 
dition for reproducing the abundance observations in CH giants near [Fe/H] = —1 (Vanture 
1992; Busso et al. 1999) and also helps in maintaining a high production of Ba over a wide 
range of metallicities. 

The trends shown in Fig. 2 can be understood through the interplay of the phenomena 
outlined above. Starting from AGB stars of nearly solar metallicity, first the s-fluence builds 
up the s-elements belonging to the Zr-peak, at neutron magic iV = 50. Then the Zr-peak 
yields decrease while the Ba-peak production increases as shown in the Figure, reaching a 
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maximum at [Fe/H] ~ —0.6. For lower metal contents, the n-flux skips Ba and feeds Pb, 
which reaches a maximum production yield at [Fe/H] = — 1 (Gallino et al. 1999). Eventually, 
also Pb decreases making the secondary nature of the s-process evident. At these very low 
metallicities, essentially all the Fe group seeds are converted to Pb. Also thanks to the 
additional effect of light (primary) neutron absorbers, below [Fe/H] ~ —2 only a very small 
s-signature is present on heavy elements. 



3. Galactic chemical evolution model 

Our model for the chemical evolution of the Galaxy is described in detail by Ferrini & 
Galli (1988), Galli & Ferrini (1989) and Ferrini et al. (1992). The model has been adopted 
to investigate different aspects of the chemical evolution of the Galaxy as well as global 
properties of external galaxies (see e.g. Molla, Ferrini, & Diaz 1997 and refererences therein). 
Here we summarize its physical characteristics relevant to our study and briefly discuss the 
modifications to the original nucleosynthetic prescriptions introduced in the present work. 

The Galaxy is divided into three zones, halo, thick disk, and thin disk, whose compo- 
sition of stars, gas (atomic and molecular) and stellar remnants is computed as function of 
time up to the present epoch i Gal = 13 Gyr. Stars are born with an initial composition equal 
to that of the gas from which they formed. The formation of the Sun takes place 4.5 Gyr 
ago, i.e. at epoch t & = 8.5 Gyr. The thin disk is divided into concentric annuli, with no 
radial flow allowed, and is formed from material infalling from the thick disk and the halo. 
In the present work we neglect any dependence on galactocentric radius in the model results 
as well as in the observational data (see Molla et al. 1997 for an application of the model 
to the study of radial abundance gradients), and we concentrate on the evolution inside the 
solar annulus, located 8.5 kpc from the Galactic center. 
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The star formation rate ip(i) is not assumed a priori, but is obtained as the outcome of 
self-regulating processes occurring in the molecular gas phase, either spontaneous or stimu- 
lated by the presence of other stars. The nucleosynthesis prescriptions are based upon the 
classic "matrix formalism" first introduced by Talbot & Arnett (1973). The mass of the 
element % ejected by a star of mass M after a time r(M) from birth is assumed to be a linear 
combination of the masses of elements j initially present in the star: 

m? Y,QuU;i'- (!) 

3 

where the matrix elements Qy depend on the mass and the metallicity of the star. 

In the framework of Galactic evolution, it is convenient to consider separately the con- 
tribution of stars in different mass ranges. The model follows the evolution of (i) single 
low- and intermediate-mass stars (0.8 M & < M < M*) ending their life as He or C-0 white 
dwarf, (ii) binary systems able to produce Type I supernovae, and (Hi) single massive stars 
(M* < M < 100 M Q ), the progenitor of Type II supernovae. The value of M+ depends 
on metallicity; we assume M± = 6 M Q for Z < 10~ 3 and = 8 M & otherwise (see for 
references Tornambe & Chieffi 1986). For each stellar group, the restitution rate W i: k(t) of 
element i is given by 

Wi, k (t)= / £Q^[M,Z(0]X^ - r{M)\^[t - T(M)]c/>(M)dM, (2) 

JM lovr ■ 

where the index k refers to the zone of the Galaxy, e.g. halo, thick disk and thin disk, and 
Xj[t — t(M)} is the abundance (by mass) of element j in the ISM at the time of birth of a 
star of mass M. The adopted initial mass function <f)(M) is discussed in Ferrini et al. (1992). 

We have updated the nucleosynthesis prescriptions of Ferrini et al. (1992) for SNII and 
SNI adopting the yields computed by Woosley & Weaver (1995) and Thielemann, Nomoto, 
& Hashimoto (1996), respectively. As an example of the model results, Fig. 3a shows the 
run of star formation rate vs. [Fe/H] in the three Galactic zones for our Standard Model 
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( corresponding to Model B of Pardi, Ferrini, & Matteucci 1995). The [Fe/H] scale here is 
indicative of a time scale, albeit a non-linear one. The corresponding values of [O/Fe] vs. 
[Fe/H] for the three zones, computed with the revised yields, are plotted in Fig. 3b, together 
with an updated compilation of observational data. We have included abundance determi- 
nations obtained from measurements of the [01] forbidden lines (Gratton & Ortolani 1986; 
Barbuy 1988; Barbuy & Erdelyi-Mendes 1989; Sneden et al. 1991; Spite & Spite 1991; Ed- 
vardsson et al. 1993) and from the OH bands in the near UV (Israelian et al. 1998). 

We see from Fig. 3a-b that the halo phase lasts approximately up to [Fe/H] <; —1.5, 
the thick-disk phase covers the interval —2.5 ^ [Fe/H] ^ —1, and the thin-disk phase starts 
at approximately [Fe/H] ^ —1.5. The slow decline of [O/Fe] for increasing [Fe/H], due to 
contribution of SNI as the evolution proceeds, is reproduced by the model over more than 
three decades in metallicity. 

Recently Israelian, Lopez, & Rebolo (1998) have measured oxygen abundances from the 
OH bands in the near UV in a sample of stars in the metallicity range —3.0 <[Fe/H]< 0.3 
(also included in Fig. 3b). They found values significantly higher than those reported by 
other authors shown in the same Figure, with a continuous increase of [O/Fe] with decreasing 
[Fe/H]. To reproduce these high values, one possibility is to increase the O contribution with 
respect to Fe of the most massive Type II SNe by a factor ~ 2. This modification, of course, 
has no consequence on the chemical evolution of the s-process elements. 



4. Results for the Galactic evolution of s- and r-elements 

In this Section we present our results for the evolution of Ba, La, Ce, Nd, Pr, Sm and 
Eu in the Galaxy, by considering separately the s- and r-contributions. Then we compute 
the abundance of these elements resulting from the sum of the two processes and we com- 



- 15 - 



pare model results with the available spectroscopic observations of field stars at different 
metallicities. 

4.1. Galactic s-process evolution 

The s-process yields discussed in Sect. 2 and shown in Fig. 2 allow us to estimate the 
chemical enrichment of the Galaxy at different times. We show in Fig. 4 the resulting Ba 
s-fraction in the thin disk, compared with spectroscopic abundances obtained from various 
observational campaigns (Gratton & Sneden 1994; Woolf, Tomkin, & Lambert 1995; Frangois 
1996; McWilliam et al. 1998; Nissen & Schuster 1997; Norris, Ryan, & Beers 1997; Jehin et 
al. 1998; Mashonkina, Gehren, & Bikmaev 1999). 

Though we have considered stars from 2 to 8 M , the mass range from which we get 
the dominant production of Ba is 2 — 4 M & (as shown in Fig. 4), well representative of the 
most common low- mass AGB stars. More massive AGBs, in the range 4 — 8 M , do not give 
a relevant contribution to Ba-peak elements (see the following discussion). 

As anticipated in Sect. 2, the direct comparisons of spectroscopic observations of AGB 
stars, solar system abundance determinations, isotopic anomalies in presolar SiC grains, with 
the predictions of stellar models constrain the mass of the 13 C pocket to within a certain 
range around a standard value. In the framework of Galactic chemical evolution, we adopt 
averaged stellar yields (shown in Fig. 1 and Fig. 2), calculated by assuming that all values of 
the 13 C pocket are equally probable within this interval. However, to test the sensitivity of 
our results to the choice of this fundamental parameter, we have also computed the Galactic 
evolution of the s-elements with stellar yields corresponding to the lower and upper bounds 
for the mass of the 13 C pocket. These aspects will be further discussed at the end of this 
subsection. 
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TABLE 2 

S-PROCESS FRACTIONAL CONTRIBUTIONS (%) AT t — t e 
WITH RESPECT TO SOLAR SYSTEM ABUNDANCES 

case T ' case 

H(b) 

case III^ KBW89^ 

Ba 10 253 80 88 

La 7 200 61 75 

Ce 8 248 75 77 

Pr 5 156 47 45 

Nd 5 178 54 46 

Sm 3 99 30 30 

Eu 0.6 20 6 3 

(a) _ 13q concentration: 0.17 times the ST case 

(b) _ 13 q concentration: 1.33 times the ST case 

(c) _ 13q pocket mass averaged 

W - Kappeler, Beer, & Wisshak (1989) 

In Table 2 we list the resulting s-fractions (with respect to the corresponding solar 
abundances) at t — t Q . In the first two columns we show the results obtained using a mass 
of the 13 C pocket scaled by 0.17 (case I) and by 1.33 (case II) from the ST case. In the 
third column we show the results for the average yields (case III). The corresponding stellar 
yields are shown in Fig. 2 (case III), and, for Ba only, in Fig. 1 (case I, II and III). We also 
compare, in the fourth column, our results with the previous estimates by Kappeler, Beer, & 
Wisshak (1989), obtained with the so called classical analysis, based on model-independent 
analytical expressions (an exponential form) for the neutron irradiation. Notice that the 
results obtained with the extreme values of the mass of 13 C are in strong disagreement 
with either the results of the classical analysis or the observed solar abundances. As an 
illustration, we see from Table 2 that at t — t Q the Ba s-fraction obtained with the average 
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yields is 80 % of the solar Ba abundance, whereas the lower and upper values of the mass of 
13 C give a Ba s-fraction of 10 % and 253 %, respectively. Notice also that the average yields 
give results in reasonable agreement with the classical analysis. 

It is evident from Fig. 4 that the s-process contribution dominates the Galactic evolu- 
tion of Ba starting from [Fe/H] ~ —1.5. At lower values of [Fe/H], independently on the 
characteristics of the chemical evolution model, the contribution of s-process nucleosynthe- 
sis rapidly decreases due to the strong dependence of stellar yields on metallicity: hence the 
contribution of low-mass AGB stars is far too low to account for the observed abundances 
in the range —3 ^ [Fe/H] <^ —1.5. In order to check the sensitivity of this result on model 
inputs, we have considered the effects of the contribution of intermediate mass stars in the 
range 4-8 M Q , adopting the yields of a 5 M© AGB model with solar metallicity (Vaglio 
et al. 1999) as representative of this mass interval. These stars activate the 22 Ne(o;,n) 25 Mg 
reaction during their TP- AGB phase more efficiently than in low-mass stars, because of the 
higher temperature reached at the bottom of the convective pulse (T max ^ 3.5 x 10 8 K). In 
contrast, the formation of a 13 C pocket is less certain, due to the reduced mass of the He 
intershell (by one order of magnitude). In this mass range the observational constraints are 
scarce and indirect: no AGB star of the Galaxy, for which we have a spectroscopic estimate of 
abundances, can be unambiguously attributed to this mass range, though in some post- AGB 
stars there is some suggestion that an intermediate mass star may have been the progenitor 
(Busso et al. 1999). Given such a lack of constraints, for these stars we have computed 
s-process nucleosynthesis induced by the 22 Ne(a,n) 25 Mg neutron source alone. We find that 
the neutron flux is too low for these stars to significantly contribute to the Galactic chemical 
evolution of the Ba-peak elements, whereas there is a consistent production of the elements 
belonging to the first s-peak: Sr, Y, and Zr. In particular, the contribution of intermediate 
mass stars only increase the Ba abundance at t — t Q by about 3%, and let the rising of 
[Ba/Fe] vs. [Fe/H] shift at slightly lower metallicity, due to the shorter stellar lifetimes. 
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The results of the Galactic chemical evolution are not particularly sensitive to the value 
of the mass contained in the 13 C pocket adopted. For example, various tests for the s-process 
nucleosynthesis have been made by progressively reducing the 13 C mass giving results very 
close to the ones illustrated in Fig. I. This is inherent to the s-process nucleosynthesis 
mechanism, which, increasing the number of available neutrons per Fe group seed (i.e. with 
decreasing metallicity), proceeds through the accumulation of neutron-magic nuclei, at N = 
50, 82, and 126, at the Zr-peak, Ba-peak and at Pb, respectively. Varying the total mass 
involved in the 13 C pocket in the present calculations by a factor of 2 is equivalent to vary by 
only 30 % the final yield, because of the overlapping mechanisms between subsequent pulses. 
This is an indication to extend the analysis to cases with somewhat higher abundances of 
13 C nuclei with respect to the upper case (1.3 x ST) considered in the present calculations, as 
suggested by the comparison of spectroscopic abundances of some s-enhanced stars (Busso 
et al. 1999). 

Despite the large number of approximations (the total amount of dredged-up material, 
which depends on the stellar evolutionary code and on the adopted law of mass loss by stellar 
winds, the parametrization on the 13 C pocket, etc.) the present results appear to reasonably 
reproduce the s-process distribution of the heavy elements, from the Ba-peak up to Pb (see 
also Gallino et al. 1998, 1999). 

4.2. Galactic r-process contribution 

As first suggested by Truran (1981), the presence of r-process elements in low metallicity 
halo stars (Gratton & Sneden 1994; McWilliam et al. 1995; McWilliam 1998) is indicative 
of a prompt enrichment of the Galaxy in these elements, possibly by early generations of 
massive stars. 
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It is possible to constrain quantitatively the enrichment of r-process elements in the 
ISM during the evolution of the Galaxy on the basis of the results presented in Sect. 4.1 
and in Table 2 for the s-process contribution at t — t & . The so-called r-process residuals are 
obtained by subtracting the s-process contribution N S /N Q from the fractional abundances 
in the solar system taken from Anders & Grevesse (1989): 



N r /N G = (N e - N s )/N e . 



(3) 



In the case of Ba we obtain a r-residual of 20%. The assumption that the r-process is of 
primary nature and originates from massive stars allows us to estimate the contribution of 
this process during the evolution of the Galaxy. In the case of Ba, for example, we have 



(4) 



Since the s-process does not contribute at low metallicity (see Sect. 4.1 and Fig. 4), for 
Population II stars we have approximately 

'Ba^ 



\OJ VO/ r , 



(5) 



which yields 

"Bal rBal TO" 

— = — + — ~log(0.2) + 0.6~ -0.1 dex, (6) 
. re J L O J l_re_ 

assuming a typical [O/Fe] ~ 0.6 dex for Population II stars. Thus, the r-process contribution 
for [Fe/H] <; —1.5 dominates over the s- contribution and roughly reproduces the observed 
values. The simple argument presented here ignores, however, the effects of the finite lifetimes 
of stars, the distribution of stellar masses, and the time dependence of the star formation 
rate. All these effects are instead taken into account in our chemical evolution model. The 
results shown in Fig. 5 for Ba and Eu confirm the estimate given above, and provide a 
strong constraint on the mass range of SNII contributing to r-process nucleosynthesis. In 
fact, the characteristic decline of [r/Fe] as function of [Fe/H] suggested by observations of 
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metal-poor stars (Cowan, Thielemann, & Truran 1991; Mathews et al. 1992; Gratton & 
Sneden 1994) can be naturally explained by a time delay between the O-rich (and partly 
Fe-rich) material ejected by the more massive SNII (M > 15 M ), and r-process material 
ejected by the lower mass SNe. This delay may reflect the actual difference in stellar lifetimes, 
but can be amplified by non- instantaneous mixing processes in the ISM. As for the scatter 
of spectroscopic data at very low metallicities, this problem will be discussed in the next 
subsection. 

In particular we have considered in Fig. 5 the following cases: (a) the full range of 
SNII, (b) SNII in the mass interval 10 < M/M Q < 12, and (c) SNII in the mass interval 
8 < M/Mq < 10. It is evident that in order to reproduce the typical increasing trend of 
[Ba/Fe] or [Eu/Fe] at [Fe/H] <^ —1.5, the production of Fe at low metallicity must have 
occurred substantially before the production of the r-process component of Ba and Eu. 
According to our model, SNII in the mass range 8 < M/M & < 10 appear to be good 
candidates for the primary production of r-nuclei, whereas a range extended to much higher 
masses seems to give results in conflict with the available observations. We note in addition 
that the identification of low-mass SNe as a possible site for the r-process is also supported 
by recent theoretical SN models by Wheeler et al. (1998), Freiburghaus et al. (1998), and 
Meyer & Brown (1997). 



4.3. s- and r-process contributions to Galactic chemical evolution 

In this Section we present our results for the Galactic chemical evolution of elements 
from Ba to Eu, based on the assumptions discussed in Sect. 4.1 — 4.2, namely that the s- 
process contribution to the production of these elements comes from 2 — 4 M Q AGB stars, 
and the r-process contribution originates from SN in the range 8 — 10 M Q . 
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Fig. 6 to 11 show the Galactic evolution of Ba, La, Ce, Nd, Pr, Sm, computed by adding 
the s- and r-process contributions (abundances are given with respect to Fe and Eu in the 
upper and lower panels respectively). Fig. 12 shows the resulting evolution of [Eu/Fe] as a 
function of [Fe/H], compared with spectroscopic abundances determined in dwarf and giant 
stars of different metallicities. 

The [element/Fe] ratios (upper panels of Fig. 6 — 11) provide information about the 
enrichment relative to Fe in the three Galactic zones, making clear that a delay in the re- 
process production with respect to Fe is needed in order to match the spectroscopic data 
in the range —3 ^ [Fe/H] <^ —2. The observations show that the [element/Fe] ratios begin 
to decline in metal-poor stars, with a sharp drop at [Fe/H] <^ —2.5. This trend can be 
naturally explained by the finite lifetimes of stars at the lower end of the adopted mass 
range: massive stars in the early times of the evolution of the Galaxy evolve quickly, ending 
as SNII producing O and Fe. Later, less massive stars explode as SNII, producing r-process 
elements and causing the sudden increase in [element/Fe]. 

Since Eu is mostly produced by r-process nucleosynthesis (94% at t — t & , see Tab. 2), 
the [element/Eu] abundance ratios (bottom panels) provide a direct way to judge the relative 
importance of the s and r channels during the evolution of the Galaxy. According to our 
model, at low metallicity the r-process contribution is dominant, and the [element/Eu] ratio 
is approximately given by the element's r-fraction computed with the r-residuals discussed in 
Sect. 4.2. On the contrary, for [Fe/H] ^ —1.5, the s-process contribution takes over, and the 
[element/Eu] ratios rapidly increase approaching the solar values. This trend is particularly 
evident in the observational data for Ba and Ce, as these elements have the largest s-process 
contribution. 

The agreement of the model results with observational data for s-process dominated 
elements (e.g. Ba and Ce) provides a strong support to the validity of the nucleosynthesis 
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prescriptions adopted in this work and AGB modelings with the FRANEC code, includ- 
ing the estimate of the s-process and 12 C-rich material dredged-up in AGB envelopes and 
eventually mixed with the ISM by stellar winds. In turn, the s-fractions obtained from our 
model at t — i , listed in Table 2, are used to normalize the r-process nucleosynthesis and 
therefore determine the element abundances at earlier Galactic epochs. In addition, a few 
points should be noticed. 

First, at high metallicities, the abundances of La, Nd, and Sm with respect to Fe are 
slightly overestimated as compared to spectroscopic observations, whereas the ratios to Eu 
are well reproduced. On the other hand, for Ba, Ce, and Pr, the predicted abundance ratios 
with respect to Fe are in good agreement with the observational data, whereas the ratios 
with respect to Eu appears to be slightly underestimated. The origin of discrepancies in the 
[element/Fe] ratio might be attributed to specific characteristics of the Galactic evolution 
model (e.g. the Fe production by SNI/SNII). The ratio [element/Eu], however, is much 
less model dependent, and discrepancies between model results and observations should be 
attributed more likely to incomplete observational data. 

Second, europium deserves a particular attention as it is one of the few r-process ele- 
ments that has clean atomic lines accessible in the visible part of the spectrum, which makes 
it an important diagnostic of the r-process history of stellar material. In Fig. 12 we show 
the observationally determined values of [Eu/Fe] in a sample of halo and disk stars. For 
[Fe/H] > —2.4, the similarity of O and Eu trends vs. [Fe/H] strongly supports the idea of 
a common production site, presumably SNII. In the same Fig. 12 we also show the results 
of our chemical evolution model, obtained following the prescriptions described in Sect. 4.2. 
The scatter in [Eu/Fe] at lower metallicities (also observed in the other Ba-peak elements) 
can be ascribed to an incomplete mixing in the Galactic halo gas, allowing the formation of 
stars super-rich in r-process elements, like CS 22892-052 (see Sneden et al. 1996). This star 
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in particular, with [Fe/H] ~ —3.1, shows r-process enhancements of 40 times the solar value 
(e.g. [Eu/Fe] ~ +1.7, and [Ba/Fe] ~ +0.9), much larger than the abundances observed 
in any normal halo star. Nevertheless its [Ba/Eu] is in good agreement with the typical 
r-process ratios we predict at this metallicity. Therefore CS 22892-052 is likely to be a star 
born in an environment more strongly polluted by supernova debris than the average halo 
gas; its peculiar abundances support the idea of an inhomogeneous composition of the halo 
at early times. 

Third, at low metallicity, the values of [element/Fe] shown in Fig.6 — 11 suggest an 
intrinsic scatter over about 2 dex. McWilliam et al. (1995) and McWilliam (1998) have 
argued that the observed dispersion in heavy element abundances at low metallicity entirely 
reflects the actual inhomogeneities in the chemical composition of the gas from which these 
extremely metal-poor halo stars were formed. Ryan, Norris, & Beers (1996) took a step 
further and speculated that the resulting chemical enrichment of the gas would be highly 
inhomogeneous and basically limited to the supernova's sphere of influence. The problem 
of the observed spread in the measured stellar abundances, in connection with mixing of 
heavy elements in the various Galactic zones is a complicated one, and in our opinion not 
fully understood. The spread observed in the Galactic disk population by Edvardsson et 
al. (1993) may result either from the multizone structure of the Galaxy (the thick and 
thin disk contribute in different ways to the local enrichment), or from a contamination of 
contiguous radial zones (see Molla et al. 1997 for details). In the halo, the non locality of 
observed populations will be again a good reason for a relevant spread in the abundance 
distribution. The question remains if these mechanisms can be responsible for the extremely 
large spread in the spectroscopic data presented by McWilliam (1998). 

We must keep in mind that models of Galactic chemical evolution often assume instan- 
taneous mixing between stellar ejecta and the ISM, and they can therefore provide only 



-24- 



spatially averaged values of element abundances as a function of time and galacto centric 
distance. While such an assumption is justified during the evolution of the Galactic disk and 
the late times of the halo phase, a different kind of approach should be taken when inves- 
tigating the earliest phases of halo enrichment, where a coupled treatment of the chemistry 
and the dynamics of the gas, both on a large scale (collapse) and on a small scale (super- 
novae), is needed for a more realistic modeling. Preliminary results concerning Ba are shown 
in Cravanzola et al. (1999), where the chemo-dynamical evolution of the Galaxy is followed 
by means of a N-body/SPH code. Similar conclusions have been advanced by Ishimaru & 
Wanajo (1999). 

On the other hand, the consistent observational scatter at very low metallicities can 
be attributed, at least in part, to uncertainties in the spectroscopic observations, which 
suffer for poor signal-to-noise ratios due to the intrinsic weakness of the sources (Gratton & 
Sneden 1994; Gratton 1995). 

The s-process contributions to the various isotopes of the elements considered in this 
paper are listed in Table 3. For example, at t = t® the model can account for 94% and 97% 
of the s-only isotopes 134 Ba and 136 Ba, respectively. The dominantly r-process Eu isotopes 
are affected by the s-process only at the 5% level. The solar distribution of s- and r-fractions 
of individual isotopes is reproduced within ~ 10%. 

5. Conclusions 

In this paper we have calculated the evolution of n-capture elements from Ba to Eu 
in the interstellar gas of the Galaxy. The input stellar yields for neutron-rich nuclei have 
been separated into their s- and r-process components. We have obtained the s-yields with 
post-process calculations based on AGB models computed with the FRANEC code. The 
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TABLE 3 

S-PROCESS FRACTIONAL CONTRIBUTIONS (%) FOR ISOTOPES FROM BA TO Eu 
AT t = t Q WITH RESPECT TO SOLAR SYSTEM ABUNDANCES 



134 Ba 


94 


141 p r 


47 


147 Sm 


20 


135 Ba 


22 


142 Nd 


93 


148 Sm 


96 


136 Ba 


97 


143 Nd 


30 


149 Sm 


12 


137 Ba 


58 


144 Nd 


48 


150 Sm 


94 


138 Ba 


84 


145 Nd 


26 


152 Sm 


21 


139 La 


61 


146 Nd 


61 


154 Sm 


0.5 


140 Ce 


81 


148 Nd 


13 


151 Eu 


6 


142 Ce 


12 


150 Nd 


0.02 


153 Eu 


5 



results of the Galactic evolution model, compared with spectroscopic observations of F and 
G dwarf stars, confirm the basic nucleosynthesis scenario outlined by Straniero et al. (1997) 
and Gallino et al. (1998). The 13 C neutron source, active during the interpulse phases of 
low-mass TP-AGB stars, accounts for most of the s-process contribution to elements from 
Ba to Eu. In the Galactic disk the abundance of Ba is dominated by the s-process for [Fe/H] 
^ —1.5. At lower metallicities, the primary r-process contribution, relatively small at the 
time of formation of the Sun (20% of solar Ba), plays a dominant role. 

Concerning the enrichment of the Galaxy in the r-process elements, the comparison 
of the predictions of the Galactic chemical model with the available spectroscopic data for 
Population II stars suggests a production from SNII slightly delayed with respect to the main 
phase of oxygen enrichment. The overall decline of [r/Fe] vs. [Fe/H] in the most metal-poor 
stars can be explained if the r-process derives from core-collapse SNe at the lowest stellar 
mass limit, around 8-10 M . 



-26- 



The Galactic evolution of the lighter s-process nuclei (e.g. Sr, Y, Zr) as well as of Pb 
will be presented in a forthcoming paper. Spectroscopic determinations of the abundance of 
these elements in Galactic halo stars are currently under way (see e.g. Sneden et al. 1998 
for recent HST data). Concerning Pb, preliminary results on the s-fraction at t — t Q and on 
its sharp increase at low metallicities, confirm the prediction by Gallino et al. (1998) that 
the production of 208 Pb (the strong component) has to be entirely attributed to s-process 
occurring in low-metallicity, low-mass AGB stars. 
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[Fe/H] 

Fig. 1. — Abundances by mass (relative to solar) of Ba in the material cumulatively mixed to 
the surface of a 2 M & star by third dredge-up episodes as function of metallicity, for different 
assumptions on the mass of the 13 C pocket. The short- and long-dashed thick lines represent 
the cases 0.17 and 1.33 times the standard value (see text) respectively, whereas the thin 
lines show intermediate cases (0.33, 0.45, 0.67, and 1 times the standard value). The thick 
continuous line represents the unweighted average of all cases shown. 
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Fig. 2. — Abundances by mass (relative to solar) of Ba, La, and Ce (top left panel), Nd and 
Pr (top right panel), Sm and Eu (bottom left panel) in the material cumulatively mixed to 
the surface of a 2 M star by third dredge-up episodes as function of metallicity. 
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Fig. 3. — Star formation rate {upper panel) and [O/Fe] {lower panel) according to our 
standard model of Galactic evolution, displayed as function of [Fe/H]. [O/Fe] predictions are 
compared with spectroscopic observations by: Gratton & Ortolani (1986) {circles); Barbuy 
(1988) {open triangles); Barbuy & Erdelyi-Mendes (1989) {four-pointed stars); Sneden et al. 
(1991) {filled triangles); Spite & Spite (1991) {filled squares); Edvardsson et al. (1993) {open 
squares); Israelian et al. (1998) {diamonds) . 
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Fig. 4. — Evolution of Ba s-fraction in the thin disk as function of [Fe/H] according to 
our standard model (solid line). Observational data are from Gratton & Sneden (1994) 
(filled squares); Woolf et al. (1995) (open squares); Frangois (1996) (pentagons); McWilliam 
et al. (1995) and McWilliam (1998) (circles); Norris et al. (1997) (filled triangles); Jehin 
et al. (1998) (filled tilted triangles; Mashonkina et al. 1999 (open triangles). Thin dotted 
lines connect stars with different abundance determinations. Errorbars are shown only when 
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reported by the authors for single objects. 
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Fig. 5. — Galactic evolution of the r-fraction of Ba (left panels) and Eu (right panels) accord- 
ing to our standard model, for different assumptions on the mass range of SNII contributing 
to the r-process nucleosynthesis (all SNII in the upper panels, 10-12 M & SNII in the mid- 
dle panels, 8-10 M SNII in the bottom panels). Lines refer to the halo (short- dashed), 
thick-disk (long-dashed) and thin-disk (solid). Data points are the same as in Fig. 4 and 12 
(errorbars are omitted for clarity). 
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Fig. 6. — Galactic evolution of [Ba/Fe] (upper panel) and [Ba/Eu] (lower panel) according 
to our standard model, including both the s- and r-process contributions. All symbols are 
the same as in in Fig. 4. Errorbars are shown only when reported by the authors for single 
objects. 
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Fig. 7. — Galactic evolution of [La/Fe] (upper panel) and [La/Eu] (lower panel) . All symbols 
are the same as in in Fig. 4. Errorbars are shown only when reported by the authors for 
single objects. 
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Fig. 8. — Galactic evolution of [Ce/Fe] (upper panel) and [Ce/Eu] (lower panel) . All symbols 
are the same as in in Fig. 4. Errorbars are shown only when reported by the authors for 
single objects. 
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Fig. 9. — Galactic evolution of [Pr/Fe] (upper panel) and [Pr/Eu] (lower panel) . All symbols 
are the same as in in Fig. 4. Errorbars are shown only when reported by the authors for 
single objects. 
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Fig. 10. — Galactic evolution of [Nd/Fe] (upper panel) and [Nd/Eu] (lower panel). All 
symbols are the same as in in Fig. 4. Errorbars are shown only when reported by the 
authors for single objects. 
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Fig. 11. — Galactic evolution of [Sm/Fe] (upper panel) and [Sm/Eu] (lower panel). All 
symbols are the same as in in Fig. 4. Errorbars are shown only when reported by the 
authors for single objects. 
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Fig. 12. — Galactic evolution of [Eu/Fe] according to our standard model, including both 
the s- and r-process contributions (even at t = teal the latter component accounts for 94% 
of the total production). All symbols are the same as in in Fig. 4. Errorbars are shown only 
when reported by the authors for single objects. 
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Fig. 13. — s-fractions at t — t e of Ba to Eu isotopes, according to our standard model. 
Values are given in percentage with respect to solar abundances. 



